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An exact theoretical framework based on Time Dependent Density Functional The- 
ory (TDDFT) is proposed in order to deal with the time-dependent quantum transport 
in fully interacting systems. We use a partition-free approach by Cini in which the 
whole system is in equilibrium before an external electric field is switched on. Our the- 
ory includes the interactions between the leads and between the leads and the device, 
ft is well suited for calculating measurable transient phenomena as well as a.c. and 
other time-dependent responses. We show that the steady-state current results from 
a dephasing mechanism provided the leads are macroscopic and the device is finite. 
In the d.c. case, we obtain a Landauer-like formula when the effective potential of 
TDDFT is uniform deep inside the electrodes. 

I. INTRODUCTION 

During the last decade, the experimental progress in the fabrication/manipulation of nano-systems 
like quantum wires and dots has enhanced the interest in the theory of quantum transport. The 
Landauer-Biittiker formalism [1,2] applies to noninteracting systems and gives the steady-state current 
of macroscopic electrodes connected via mesoscopic (or nanoscopic) constrictions. Techniques based 
on non-equilibrium Green functions [3,4] provide a natural framework for obtaining the Landauer 
formula starting from a microscopic theory. The first attempt was made by Caroli et al. [5,6] who 
consider the two leads as isolated subsystems with different chemical potentials in the remote past. 
The current will flow through the system once the contacts between the device and the leads have 
been established, see Fig. 1(a). This approach is based on a fictitious partitioning since in a real 
experiment the whole system is in thermodynamic equilibrium before an external bias is applied deep 
inside the electrodes. Later on this partitioned scheme was adopted by Meir and Wingreen [7] to 
obtain the steady-state current through an interacting central device. Shortly after, Wingreen et al. 
[8-10] generalised the Meir- Wingreen formula to time-dependent phenomena. 

Despite the important results mentioned above, the partitioned scheme has several drawbacks. 
First, it is difficult to include the interactions between the leads and between the leads and the device. 
These interactions are responsible for the establishing of dipole layers and charge transfers which 
shape the potential landscape in the device region and influence the contacts. Second, there is a 
crucial assumption of equivalence between the long-time behaviour of the 1) initially partitioned and 
biased system once the coupling between the subsystems is switched on and 2) the whole partition-free 
system once the bias is applied. Third, there is no well defined prescription which fixes the initial 
equilibrium distribution of the isolated central device; this makes the transient current difficult to 
interpret. 

A conceptually different approach from the one by Caroli et al. has been introduced by Cini [11], 
see Fig. 1(b). Here, the system is contacted and in thermodynamic equilibrium before an external 
time-dependent disturbance is switched on. Cini developed the general theory for the case of free 
electrons with his focus on semiconductor junction devices. In this letter we extend the Cini theory to 
fully interacting (electrodes and device) systems. Our exact theoretical framework is based on Time 
Dependent Density Functional Theory (TDDFT) and is well suited to describe the steady state as 
well as a.c. and other time- dependent current responses. 
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FIG. 1. (a) Schematic representation of the partitioned approach by Caroh et al. On the left the electrodes are 
disconnected and in thermodynamic equilibrium at two different chemical jjotential. On the right the contacts are 
established and electrons start to flow, (b) Schematic representation of the partition-free approach by Cini. The whole 
system is in thermodynamic equilibrium and it is characterised by a unique chemical potential. Then, a bias is applied 
deep inside the electrodes. 



II. GENERAL FORMULATION 



To be specific let us focus on the Coulomb interaction and on paramagnetic systems. In the case we 
only ask for the time-dependent density n(r,t) the original Density Functional Theory [12,13] and its 
finite-temperature generalisation [14] has been extended to time-dependent phenomena [15,16]. This 
theory applies only to those cases where the external disturbance is local in space, / dr U {r , t)n{r , t) . 
Runge and Gross [15] have shown that one can compute n{r,t) in a one-particle manner using an 
effective potential 

[/^ff(r,i) = f/c(r,i) + ^;xc(r,t). 

Here, Uc = U+Vu + Vnis the classical electrostatic potential and it is given by the sum of the external 
potential U, the Hartrcc potential Vh and the Coulomb potential from the nuclei Vn- Furthermore. Uxc 
accounts for exchange and correlations and is obtained from an exchange-correlation action functional, 
^^xc(r, = dA:^c['>T]/5n{r,t) (as pointed out in Ref. [17], the causality and symmetry properties require 
that the action functional Axc[fi] is defined on the Keldysh contour). Thus, the effective one-body 

Hamiltonian reads = J dr ^— + ?7'^*(r, i) . Here and in what follows, we use boldface 

for operators in one-particle Hilbert space, and ?i is 1 in our units. 

It is convenient to define the projectors Pa = dr\r){r\ onto the left or right leads (a = L,R) or 
the central device (a = D) and to introduce the notation Oa/) = PaOPjs, where O is an arbitrary 
operator in one-body space [18]. The one-particle scheme of Time Dependent Density Functional 
Theory (TDDFT) corresponds to a fictitious Green function Q{t;t') which satisfies a one-particle 
equation of motion 



(1) 
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In eq. (1) £ = T-C^l + '^'dd + ^/ffl is the uncontacted one-body Hamiltonian of TDDFT while 
V = 'W^ — £ accounts for the contacting part. The fictitious Q wiU not in general give correct 
one-particle properties. However, by definition, gives the correct density 

n{r,t) = -2i{r\g<{t;t)\r) 

(where the factor of 2 comes from spin). Also total currents are correctly given by TDDFT. If for 
instance Ja is the total current from a particular region a we have 

Mt) = -e j dr jn{v,t), (2) 

where the space integral extends over the region a (e is the electron charge). The TDDFT gives the 
exact density, and from eq. (2) we see that it also gives the exact longitudinal current density and 
total current from subrcgions. 

Since V lr = V rl — 0, from eq. (1) and eq. (2) the current from the a = L,R electrode to the 
central region can be written as 

Ja{t) = 2e^[tv{g<{t;t)Vc,D}] , a = L,R. 

Without loss of generality we may assume that the partition- free system is in equilibrium when t <0, 
i.e., £(t < 0) = £'^ . For the noninteracting system of TDDFT everything is known once we know 
how to propagate the one-electron orbitals in time and how they are populated before the system 
is perturbed. The time evolution is fully described by the retarded or advanced Green functions 
Q^'^, and the initial population at zero time, i.e., by 5^(0,0) = if{£° + V), where / is the Fermi 
distribution function [since £^ + V is a matrix, so is f{£ + V)]. Then, for any t,t' > we have 
[19,11,20] 

g<{t, t') = g^{t, o)a< (0, o)a^(o, t') = i o)/(£° + F)e^(o, f), (s) 

and the total time-dependent current can be rewritten as 

Mt) = -2e 9 [tr {g«(t, 0)/(£° + y)g^(0, t)VaD}] • (4) 

The above equation is an exact result. For noninteracting electrons, eq. (4) agrees with the formula 
obtained by Cini [11] in the partition- free scheme. Indeed, the derivation by Cini does not depend on 
the details of the noninteracting system and therefore it is also correct for the Kohn-Sham system, 
which however has the extra merit of reproducing the exact density. The advantage of this approach is 
that the interaction in the leads and in the conductor are treated on the same footing via self-consistent 
calculations on the current-carrying system. It also allows for detailed studies of how the contacts 
influence the conductance properties. We note in passing that eq. (4) is also gauge invariant since it 
does not change under an overall time-dependent shift of the external potential which is constant in 
space. Neither, is it modified by a simultaneous shift of the classical electrostatic potential and /j. for 
t < 0. 



III. ASYMPTOTIC CURRENT 



Once the effective potential of TDDFT is known, the problem is reduced to a fictitious one-particle 
problem. For any such one-particle problem we now show that the partitioned scheme by Caroli et 
al. and the partition-free scheme are equivalent in the long-time limit. To this end, we consider the 
Dyson equation that relates a noninteracting system contacted all the time (Green function g) to the 
corresponding uncontacted one (Green function g). Both systems are exposed to the same, possibly 
time-dependent, potential for i > 0. As is well known, the noninteracting Green function g projected 
in a subregion a = L, R oi D can be described in terms of self energies which account for the hopping 
in and out of the subregion in question. Considering the device Green function, the self energy can 
be written = J2a=L,R ' X = R, A, < or >. Here, is given by 

i:^{t;t') = VDa9'^{t;t')V^D, 



3 



where g satisfies eq. (1) with V = 0. The g can be expressed in terms of the one-body evolution 
operator S{t) which fulfils iS{t) = £{t)S{t) and <S(0) = 1. The retarded and advanced components 
are Since also the uncontacted system is initially in equilibrium 

wehaveg<{t;t') = ig^{t;0)f{£°)g^{0;t') [cf. eq. (3)]. 
The two schemes are equivalent if 

lim t') = hm S^(t'; t) = (5) 

t — ^oo t — ^oo 

for any non-singular V. Indeed, from the equation of motion (1) and eq. (5) one can verify that the 
projected Green functions satisiy 

lim g%{t; t') = lim g^^{t'; t)Vo.D = 0. (6) 

t — ^CxD t — ^OO 

Making use of the above relations one finds [20] that the asymptotic {t oo) time-dependent current 
of eq. (4) becomes 



Mt) = -2e3 [tr{g^(t,0)/(£°)g^(0,t)V„D}] . (7) 



Thus, the long-time limit washes out the initial effect induced by the conducting term V . Moreover, 
due to eq. (6), the asymptotic current is independent of the initial equilibrium distribution of the 
central device. Expressing and Q^^ in terms of and Q^d respectively, eq. (7) can be 
rewritten in a more familiar form 

Mt) = 2eJ dm [tr{ggc(<,t)S<(?,t) +g<c(t,t)S^(f,t)}] , (8) 
where the asymptotic relation (t, t' — > 

oo) QoDit^t') = J dtdt' g^ij{t,t)i:<{tJ')g^jj(t',t') has been 
used. Eq. (8) is valid for interacting devices connected to interacting electrodes. It provides a useful 
framework for studying the transport in interacting systems from first principles. Our results can be 
applied both to the case of a constant (d.c.) bias as well as to the case of a time-dependent {e.g., a.c.) 
one. 

To summarise, if the retarded/advanced self energy satisfies eq. (5), the exact formula in eq. (4) 

reduces to eq. (8) in the long-time limit. For noninteracting electrons the Green function g of TDDFT 
coincides with the Green function of the real system and eq. (8) agrees with the formula by Wingreen 
et al. [8,9] 

IV. STEADY STATE 

Next, we wish to investigate the steady state, i.e., limt^oo£{t) = £°° = const. In this case it 
must exist a unitary operator S such that limt^oo = e~^^ * S. Then, in terms of diagonalising 
one-body states \(l)ma) ^'aa ^i^h eigenvalues e^^ we have 

^<{t;t') = iJ2 e-*"<'*-^"'«*V^„|Ca>(Cal/(^°)IC'a)(C'J^ai^, 

m.m' 

where £° = 5 £'''3^ . For t,t' oo, the left and right contraction with a nonsingular V causes a 
perfect destructive interference for states with \e^^ — e^,^\ > l/{t + t') and hence the restoration of 
translational invariance in time 

(^) t ) = i ^ ^ fma^moi^ mai )^ ^9) 



where = (0mal/(^°)l<^ma) while Tma = VDa\(l)ma){<l^ma\yaD [21]. The abovc dephttsiug mech- 
anism is the key ingredient in the developing of a steady state. Substituting eq. (9) into eq. (8) we 
get the steady state current 

J„ = -2eJ2fmf3 Wf3tT{T^rno.^[g%(e^J]\+tT{gl^ie:^0)TrnfsGU^^^ 

m/3 ^ ^ ^ ^ (10) 
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with gS'^(e) = [e- - S^'^(e)]-^ Using the equahties ^[GId\ = [QId " QdoV^^ and [^Bd - 
Qdd] = [Qdd - Qdd] together with 

ma 

and 9[S^(e)] = TrJ^m^i^ ~ ^ma)^rna, the steady-state current in eq. (10) can be rewritten in a 
Landauer-hke form 

Jr = —^''^[fmL%nL — fmR%nR] = —Jl- (H) 

m 

In the above formula 7^_r = J2n ^mR and %nL = X)„ are the TDDFT transmission coefficients 
expressed in terms of the quantities 

"^"a = 27r(5(e5^jj — e^)tr |5^j3(e5^Q,)r„iQ,^^£,(eJ5^)r„^| = 7^". 

Despite the formal analogy with the Landauer formula, eq. (11) contains an important conceptual 
difference since is not simply given by the Fermi distribution function. For example, if the induced 
change in effective potential varies widely in space deep inside the electrodes, the band structure 
may be completely different from that of However, if we asymptotically have equilibrium far 

away from the central region, as we would expect for leads with a macroscopic cross section, the 
change in effective potential must be uniform. To leading order in 1/N we then have 

Saa{t)=Sl^ + i^Uf{t), (12) 

and 5^ = + Ai7°^. Hence, except for corrections which are of lower order with respect to the 
system size, = and 

fma = fi^ma ~ ^^a,oo)- 

We emphasise that the steady-state current in eq. (11) comes out from a pure dephasing mechanism 
in the fictitious noninteracting problem. The damping effects from scatterings is described by A^c and 
Uxc- Furthermore, if eq. (12) holds, the current depends only on the asymptotic value of the effective 
potential, U^^{r, t — > oo). However, U'^^{r, t — > oo) may depend on the history of the external bias U. 
Thus, the steady state current of fully interacting systems may be history dependent. In the case of 
Time Dependent Local Density Approximation, the exchange-correlation potential Vy^c depends only 
locally on the instantaneous density and has no memory at all. If the density tends to a constant, so 
docs the effective potential ^ which again implies that the density tends to a constant. Owing to 
the non-linearity of the problem there might still be more than one steady-state solution or none at 
all. 

The Landauer formula in eq. (11) corresponds closely to the result obtained by Lang and coworkers 
[22,23]. In their approach, the continuum is split into left- and right-going parts which are populated 
according to two different chemical potentials. The density is then calculated self-consistently. In our 
language this corresponds to writing the uncontacted g^(e) as 

g<{e) = 2mJ2Ue^JS{e - e~J|Ca)(Cal 

ma 

in terms of Fermi functions fa with chemical potential yU^ = /U -|- AC/^*^. The chemical potentials 
for the two leads differ, and the final result is independent of the chosen chemical potential for the 

device due to eq. (5). When we apply 1 -|- Q^V = Q^[g^]~^ to an unperturbed state \4>ma)^ it is 
transformed to an interacting, i.e., contacted eigenstate IV'ma}- Above the conductance threshold, 
states originating from the left lead become right-going scattering states, and states from the right 
lead become left-going scattering states. In addition, fully refiected waves and discrete states may 
arise which contribute to the density but not to the current. Thus, 

g<{e) = [1 + g%)V] g<{s)[l + Vg^ie)] = 27r» ^ /„(e-J5(£ - e^aM^^JO^mJ- 

ma 
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Lang et al. further approximate exchange and correlation by the LDA and the leads by homogeneous 
jellia, but apart from these approximations it is clear that their method implements TDDFT, as 
described above, in the steady state. It is also clear that the correctness of the Lang's approach relies 
on eq. (5) and eq. (12), as we have seen before. The equivalence between the scattering state approach 
by Lang et al. and the partitioned non-equilibrium approach used by Taylor et al. [24,25] has also 
been shown by Brandbyge et al. [26] 

V. CONCLUSIONS 

In conclusion, we have proposed a partition-free scheme based on TDDFT in order to treat the 
time-dependent current response of fully interacting systems. Our formally exact theory is more akin 
to the way the experiments are carried out and allows us to calculate the physical {i.e., measurable) 
dynamical current responses. Among the advantages we stress the possibility to include the interac- 
tions between the leads and between the leads and the device in a very natural way. We have shown 
that the steady state develops due to a dephasing mechanism without any reference to many-body 
damping and interactions. The damping mechanism (due to the electron-electron scatterings) of the 
real problem is described by v^c- The steady-state current depends on the history only through the 
asymptotic shape of the effective TDDFT potential C/°^ provided the bias-induced change AU^^^^ is 
uniform deep inside the electrodes. (This is the anticipated behaviour for macroscopic electrodes.) In 
such systems, we have shown that the nonlinear steady-state current can be expressed in a Landauer- 
like formula in terms of fictitious transmission coefficients and one-particle energy eigenvalues. Our 
scheme is equally applicable to time-dependent responses. Clearly, its usefulness depends on the qual- 
ity of the approximate TDDFT functionals being used. We are presently investigating how different 
approximations influence the current response in model systems. 
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